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We show that the isocurvature perturbations imply that the gravitationally produced superheavy 
dark matter must have masses larger than few times the Hubble expansion rate at the end of 
inflation. This together with the bound on tensor to scalar contribution to the CMB induces a 
lower bound on the reheating temperature for superheavy dark matter to be about 10 7 GeV. Hence, 
if the superheavy dark matter scenario is embedded in supergravity models with gravity mediated 
SUSY breaking, the gravitino bound will squeeze this scenario. Furthermore, the CMB constraint 
strengthens the statement that gravitationally produced superheavy dark matter scenario prefers a 
relatively large tensor mode amplitude if the reheating temperature must be less than 10 9 GeV. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The primordial density perturbations which provide initial conditions for structure formation Q can be classified as 
two types: curvature and isocurvature. Isocurvature perturbations are defined to be fluctuations in the composition 
of the energy density. Hence, the curvature remains fixed while the relative contributions of various fluid elements 
composing the total energy density changes. More explicitly, the isocurvature perturbation due to a cold dark matter 
particle species X can be written as 

Sn x Sn~ 

S s = j 1 

nx n 7 

where rt 7 is the density of photons and nx is the density of the dark matter particles. As can be seen in the appendix, 
generically, isocurvature perturbations are almost always synonymous with entropy perturbations. Although the 
isocurvature perturbations do not give rise to potential energy during radiation domination, after matter domination 
occurs, the isocurvature perturbations also become source for gravitational potential energy leaving an imprint on the 
CMB. 

In the late 80's and early 90's, isocurvature and curvature perturbation models were seen to be competing models 
for the theory of initial conditions for structure formation. Today, from the precise measurements of the CMB acoustic 
peaks and initial measurements of the TE correlations 0, Q , most people agree that the dominant contribution to 
the primordial density fluctuations arise from curvature perturbations. Nonetheless, as demonstrated by many recent 
works [3, la, 111 III IE 111 lill j an order 10% contribution to the CMB power spectrum from isocurvature is possible. 
Perhaps more importantly, we generically expect from theoretical considerations a nonvanishing contribution to the 
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CMB from isocurvature perturbations. (In other words, whenever there is more than one energy component in the 
universe, there is always some amount of isocurvature perturbations.) 

In the most popular WIMP CDM (cold dark matter) picture in inflationary cosmology, the isocurvature pertur- 
bations are expected to be small because CDM abundance today is determined from a freeze out process starting 
from chemical equilibrium initial conditions. Since the chemical equilibrium is with radiation, which, by definition, 
has overdensitics defined by curvature fluctuations during radiation domination, CDM in this case will naturally have 
suppressed isocurvature component since 

Snx Sn~ ( 3ST 
nx n-y 1 

where X denotes the dark matter density, 7 denotes radiation, and T denotes the equilibrium temperature. 

However, for superheavy dark matter which never equilibrates, the isocurvature perturbations are generically not 
negligible. Furthermore, because, for gravitational particle production scenario, the superheavy dark matter ho- 
mogeneous density is determined by the same inflationary dynamics as the isocurvature perturbations themselves, 
one obtains a correlated constraint for the superheavy dark matter abundance and the isocurvature perturbation 
amplitude. In this paper, we compute the generic feature of this correlated constraint. 

The gravitational production of superheavy dark matter has been well studied in the past 0, 0, 0, 0, 0] . In 
this scenario, the nonadiabatic change in the particle dispersion relationship due to a relatively sudden transition 
from a quasi-de Sitter phase of inflation to a small power law expansion phase of FRW cosmology causes particle 
production. The only non-generic feature of this dark matter scenario is the existence of stable particles with heavy 
mass. For 10 11 — I0 13 GeV particles, up to dimension 8 decay operators needs to be absent to have the particle be long 
lived enough to be CDM. Nonetheless, there are many par ticle physics and string theory motivated candidates for 
superheavy dark matter (see for example |16l Il7l Il8l Il9l l2Cl 12 lll22L 123. l24l l25l |26|). Hence, what makes this scenario 
exciting is that if we can measure the signatures of superheavy dark matter, we would have a probe of the very early 
universe where the physics far beyond the standard model of particle physics was important. (Indeed, it may even 
probe string theory.) This paper serves as one step towards identifying the signature of superheavy dark matter on 
the CMB. 

In particular, we investigate how much superheavy dark matter isocurvature perturbations are produced at around 
50 efolds before the end of inflation and how those isocurvature perturbations will show up on the CMB. We find that 
the mass of the superheavy dark matter cannot be too small if one is to avoid overproducing isocurvature perturbations 
inconsistent with the CMB data. More specifically, the constraint is 

^>0(5), (3) 

where mx is the mass of the superheavy dark matter and H e is the Hubble expansion rate at the end of inflation. 
One implication of this is to strengthen the preference of large tensor perturbations for gravitationally produced 
superheavy dark matter scenario. More specifically, if the reheating temperature is less than 10 9 GeV, then the tensor 
to scalar power ratio is bounded from below by 

^ > ICr 3 . (4) 
Pr ~ K ' 

Another implication of Eq. |(3J is that there is a lower bound on the reheating temperature for superheavy dark mat- 
ter scenario to be viable since the number density of particles produced falls off as mx / H e and increasing the reheating 
temperature increases the final number density. Putting this together with the bound on the tensor perturbations 
from inflation, we find the approximate requirement 

Trh > 10 7 (^\ GeV (5) 

in order to have the superheavy dark matter be the CDM while not contributing too much to isocurvature pertur- 
bations where rs represents the bound on the ratio of tensor to scalar perturbation power. Hence, we have a severe 
restriction on this scenario coming from the gravitino bound (see for example |27l l28l l29l l3Ct l3lL 1321 133L l34 | 1 if this 
scenario is embedded in gravity mediate SUSY breaking scenario (for an introductory review of status of low energy 
supersymmetry, see for example, |35jV Note that given that both the superheavy dark matter and the gravity waves 
are produced by inflationary energy density, it is natural that there is a bound of the form Eq. (JjJJ. What is not 
obvious about Eq. J5J is how the bound on isocurvature perturbations make Eq. J5J independent of the mass of the 
superheavy dark matter. The explanation of this is one of the main results of this paper. 
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In the present work, we will consider superheavy dark matter model of massive scalar field X minimally coupled to 
gravity without any other interaction for X. The results are not likely to change drastically with fcrmionic superheavy 
dark matter since none of the physics relies upon the existence of a vacuum expectation value of the dark matter field 
nor spin statistics. Furthermore, we restrict ourselves to minimal coupling in this paper. For any fixed dark matter 
mass, nonminimal coupling will change the total amount of dark matter produced, the spectral distribution of the 
dark matter, and the dark matter's stress tensor correlator. The most important effect is likely to come from the dark 
matter's stress tensor correlator. Because a coupling of the form £,cf> 2 R (where R is the Ricci scalar) effectively shifts 
the mass of the superheavy dark matter as m 2 x — » m 2 x + 12£H 2 (where H is the expansion rate during inflation), a 
positive £ results will relax the mass bound shown in Eq. © , possibly to an extent where there will be no constraint 
from isocurvature perturbations. On the other hand, a negative value for the nonminimal coupling will make the 
bound more stringent. 

The order of presentation is as follows. In Section 2, we make general estimates to establish the results associated 
with Eqs. (0, (@J, and (JSJ). Following that, we review the basic physics of CMB and its relationship to isocurvature 
perturbations. We then compute the isocurvature power spectrum (essentially stress tensor correlation function) in 
Section 4. In Section 5, we present numerical results with a special emphasis on a particular inflationary model 
as a check of our general argument. Finally, we conclude in Section 6. In the appendix A, we give an elementary 
thermodynamic discussion of isocurvature perturbations to remind the reader of its physics. The technical details of 
the particle production computation for the toy model used in Section 5 is presented in Appendix B. 



II. GENERAL ARGUMENTS 

It is possible to semi-quantitatively describe the mass bound coming from isocurvature perturbations as follows. 
Field fluctuations are typically characterized by its iV-point moments, and the most elementary nontrivial moment 
for perturbative theories is typically the two point function. The power spectrum of any two point function in a 
translationally invariant space is defined as 

Px{k) = &J d3 ^*W)*(y)), (6) 

where r is defined to be r = x — y. As we will discuss below, the isocurvature perturbations that we are concerned 
with are characterized by the power spectrum 

P Sx (k) = ^J d 3 re-^{6 x (x)6 x (y}), (7) 

where Sx = 6px / Px is the CDM (superheavy dark matter in our case) energy overdensity. 

As an order of magnitude estimate, we can represent the energy overdensity field as Spx ~ m 2 x X 2 , which gives 

Pi* ~ / e-^(X 2 {x)X 2 {y)) ~ / d\ e^? (X{x)X{y)) 2 . (8) 

Px J Px J 

Hence, in the long wavelength limit, we can approximate this as 



Px 



For scalar fields X in FRW spacetimes in general, if V"(X) 3> H 2 (where V is the scalar field potential and H is 
the expansion rate), we know that the correlation function will behave as in flat Minkowski space. In that case, we 
know from Minkowski space field theory that the two point correlation functions receive most of its power from short 
distances. We also know that for long wavelengths (k — > 0) in de Sitter (dS) space, when V" <C H 2 , that 



H\ 2 ( k 



p x~[r«){is)- < 10) 

where (k/aHy represents the factor that is breaking scale invariance due to the fact that V"(X) ^ 0. 1 Hence, using 
the Minkowski intuition and the fact that there is generically large power contribution in the infrared for V" < (due 



1 It is important to emphasize that the form of Eq. 1101 is only valid when V" /H 2 <C 9/4 in the minimal coupling scenario. On the other 
hand, for nonminimal coupling, this form is valid only if V" /H 2 9/4 — 12£. 
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to instability) , we can estimate the power index of Eq. i|10|) as 



3 = Q 



H 2 ' 



where Q > is a positive order unity dimensionless number. We thus find 

,4 / tt\ 4 / , x 2QV"/H 2 



P Sx (k) 



(ii) 



(12) 



Since inflation eventually ends and since the gravitational particle production does not really "occur" until the 
end of inflation, we have to assign the time at which we evaluate the time dependent quantities carefully. Since the 
isocurvature fluctuations are generated as the scales leave the horizon early in the inflationary epoch, we assign H to 
the expansion rate at about 50 efoldings before the end of inflation. On the other hand, the scaling of the wave vector 
continues until the end of inflation at which time particles are nonadiabatically produced. After the nonadiabatic 
particle production occurs, the power spectrum should remain approximately time independent in the long wavelength 
limit owing to the arguments of causality and adiabatic evolution. Hence, we write 



P Sx (k) 



x 



Hi 

2vr 



a e Hi 



2QV"/H] 



(13) 



where Hi represents the expansion rate at about 50 efolds before the end of inflation and p e represents the CDM 
energy density at the end of inflation. In Section 4 we will compute this more carefully and find that although this 
gives the right order of magnitude, there is comparable contribution coming from the kinetic terms neglected in this 
estimate. 

The dark matter energy that appears in the denominator of Eq. ()13fl is not easy to estimate for minimal coupling 
because significant particle production contribution comes from the infrared which have ambiguous vacuum boundary 
conditions. As shown in the appendix, neglecting this infrared contribution, one finds 



0.7 x 10 _3 m Wm^m*) 3 / 2 exp [ 5mx ) < Px (t e ) < 1.6 x lO^mx^mx) 372 exp [ 4 " 2mx ) 



for V = Tjirficf) 2 potential. More generally, we expect that 



p e ~ 10 3 ?n x (ff e mx) 3/2 exp 



-2m x \ 
H e J 



Hence, we generically find 



PS X 600 to X \H e ) \d e Hi 



2Qm x /H, 



6XP \-H-) ' 



(14) 



(15) 



(16) 



where we have assumed V" w m x consistently with the particle production. Although naively, one would guess that 
the exponential would dominated on the right hand side of Eq. (|16|l , in reality the fc-dependent factor dominates since 
for large scales of interest for CMB, we have 



In 



k 



a e Hi 



In 



ap H 
a e Hi 



-55 + In 



\ 10 9 GeV 



-1/3 



V K 



10 15 GeV 



-1/6 



which implies that the fc-dependent term dominates: 

Hi 



Ps x (a H ) ~ 600- 



mx 



HA' 

He) 



exp 



f 4mx 

V He 



110Q 



(17) 



(18) 



Hence, we see the inflationary model dependence of the power sepctrum clearly in the ratio of Hi/H e (ratio of 
the expansion rates at the beginning and end of inflation) and mx/Hi. Despite this model dependence, we would 
generically conclude mx / Hi > O(l). More precisely, since the isocurvature contribution must not exceed the adiabatic 
contribution to Ci, we have for small I ~ 2 



C, 



(X) 
1=2 



0.1P Sx (a Q H ) <10 



-10 



(19) 
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which gives a bound of 



H e ~ VQx\l 1+ Qx 2 + Qx 2 ' (20) 



where x = 7H e /Hi is close to unity. If we for example take Q = 2/3 and x = 1, we obtain mx/H e > 6. 

Since this is similar to the result of our more careful computation, we will use this bound as our example for the 
rest of this section. However, more generally, the bound is m x /H e > 0(5). 

Now, according to our estimate of particle production shown in the appendix, the dark matter density today is 
given by 

W « 4.31 x 10 -5Zi^^4 „ i 3 f^^YJ^ ( ?*) e -2".WH, (21) 
T pj(ie) Vl0 13 GeV/ 10 9 GeVV^e/ 

This implies that if mx / H e > 6, we must increase the relic abundance by increasing mx from 3 x 10 13 GeV while 
keeping the exponential fixed or increase Trh or increase both mx and Trjj. Note that H e dependence only comes 
through mx/H e . 

Keeping Trh fixed, increasing mx from 10 13 GeV for a fixed ratio mx/H e > 6 is restricted by a bound on 
the amplitude of the gravitational wave (limiting H e ). The CMB temperature fluctuation normalization Q fixes 
Piz w 2 x 10~ 9 and since the tensor perturbation power spectrum is 

ft -(=aS))'- (22) 

where Mp = 2 x 10 18 GeV is the reduced Planck mass, we have the tensor to scalar ratio 



which may be a detectable amplitude for Hj = 10 14 GeV |3a, l37| . If the experimental bound on the tensor to scalar 
ratio is r$ {i.e., Pt/Pr. < r s)? we have 

Hi < 7 x 10 14 Vr^GeV, (24) 

and consequently, we cannot raise H e much above 10 14 GeV. 2 Therefore, by Eq. I|23|) and recalling that we are for 
the moment focusing on fixing mx/H e > 6 such that flxh 2 ~ 0.1, we find the bound on the superheavy dark matter 
mass of 



6<^<12 



1 i ( He \ , 1 i Trh 



1 + 777 ln -TF? r S + 777 In 



24 \Hj J 24 10 9 GeV 



(25) 



where we have approximately solved a transcendental equation by iteration. This corresponds to an upper bound on 
the dark matter mass of 



m x < 10 



15 



* rs T «» 



24 \Hj J 24 10 9 GeV 



(26) 



Now, suppose we want to find the minimum allowed Tru while maintaining m x /H e > 6. We can then solve for 
T RH in Eq. J2U with fl x h 2 = 0.1 and use Eq. to obtain 

,0.2 / Hj N 2 



TRH ^ w ^[rt) GeV < (27) 

where the 0.2 factor would be much smaller if we were not restricted to mx / He ^ 6. This may be dangerous for the 
gravitino bound if the scenario is embedded in gravity mediated SUSY breaking scenario. In our explicit computations 
for V((p) — ^m^cf) 2 model of slow roll inflation, we obtain a bound of > 10 8 GeV. 



2 Eq. 1241 with rg = 1 is about the same bound one obtains in single field slow roll inflationary scenarios. 
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From Eqs. (126[) and 1|27|) . we see that if we tighten the bound on the tensor amplitudes (by making rs smaller), the 
superheavy dark matter scenario can be severely restricted. In particular, the linear dependence of Eq. i|27|) on l/r$ 
can change the reheating lower bound significantly. 

Note if we require that Trh < 10 9 GeV (for example, because we want to evade the gravitino bound), we have from 
Eq. J2U and Eq. JUJ, that 

Since Vixh 2 > 0.1, we conclude mx >3x 10 13 GeV, which implies 

H e > 3 x 10 12 GeV, (29) 

or according to Eq. (|23|1 that 

?L> W -3(lL) 2 . (30) 

Contrast this with the case in which we do not impose mx / H e > 6. In that case, we can fix the value of the relic density, 
Eq. H21f) . to be flxh 2 =0.1, impose Trh < 10 9 GeV, and minimize the value of H e to obtain H e > 3 x 10 n GeV, 
where the minimzation of H e occurs for mx/H e = 5/4. Since H e is one order of magnitude smaller, the gravitational 
wave signal in this case is bounded from below by a number that is 100 times smaller than Eq. (|30(l . Hence, the 
isocurvature perturbation constraint strengthens the preference of large tensor perturbations for the gravitationally 
produced superheavy dark matter scenario. 

III. CMB TEMPERATURE FLUCTUATIONS AND ISOCURVATURE PERTURBATIONS 

A. Dynamics in the tight-coupling regime 



Here we follow the treatment of the Mukhanov ,38j (for other good analytic treatments, see |3Jj|40j). The back- 
ground FRW metric is taken with flat spatial sections and conformal time, whose time derivative is denoted with a 
prime ('). The scalar metric perturbation parameterization is chosen to be 

ds 2 = a 2 [(1 + 2$)dr 2 - (1 - 2<$>)5 tk dx t dx k ] . (31) 

Assuming that the cold dark matter (denoted with subscript CDM) is decoupled from the baryon-photon plasma, 
the dark-matter stress-energy conservation equations lead to 

(6 X - 3$)' + auh DMli = (32) 
[a(8 x - 3$)']' - aA$ = 0. (33) 

where u x is the velocity of the species X. 

With the approximation of nonrelativistic baryons, the conservation equation Tq. q = for baryons leads to 

(* 6 - 3$)' + auii = 0, (34) 
which corresponds to the conservation of baryon number. The photon energy conservation leads to 



(S 1 - 4$)' + -au^t = 0. (35) 



Assuming that the photons and baryons are tightly coupled, we set u 1 = u\ = uL, which from Eqs. (|34|) and l|35|) 
leads to 

3 

-<5 7 — 5b = constant. (36) 
With this tight-coupling approximation, the T". a — of the baryon-photon stress tensor leads to 

\ [a 5 ( Pb + Pl + P 1 + P b )u\] ' - Jr/Au', + A5(P 7 + P b ) + (e 7 + e b + P 7 + P b )A$ = 0, (37) 
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where r\ is the viscosity: 



V = y^P7 T 7> ( 38 ) 



with r 7 the mean free path of the photons. Assuming no isocurvature perturbations from the baryons, we use 

S b = ^<5 7 (39) 



and arrive at 



with the sound speed 



= , o? /,. ■ ( 41 ) 



3 1 + 3p fe /4p, 

Finally, the 00 component of the Einstein's equations is 



2 

a _ , _ / a 



A*-3-$'-3(-) d>=i^a 2 ^ x + ^). (42) 



a \ a 



A/ 2 A 3d 



The fields that need to be determined are {Sx,5 7 ,§} and the three independent equations are Eqs. (|33(l . (|40|) . and 
1)420. (The only remaining equation simply determines u l CDM .) 

As with any differential equations, these need boundary conditions. The CDM fluctuation initial spectrum is 
determined by a quantum computation of (SxSx) which is assumed to evolve collisionlessly at least until decoupling. 
Similarly, quantum computation of adiabatic perturbations provide another set of boundary conditions. Finally, 
since the 00 component of Einstein's equation, Eq. (|42|l . is a parabolic equation for $, we only require one boundary 
condition for <f>. The fact there is a growing solution supported by a source obviates the need for a boundary condition 
for $ for the cases of our interest. 



B. Relationship between CMB temperature and isocurvature fluctuation 

Let the temperature of the CMB photons be a field T + AT(r, x,p) where p = p l c /\p c \ corresponds to the direction 
of photon propagation (the subscript refers to coordinate momentum). We would like to find AT(ro,x,p), the 
temperature fluctuation of photons today, given AT(ri,x, p), the temperature fluctuations of photons at the time of 
last-scattering surface. We can define the relativistic phase space volume as 

\g a b - n a n b \d 3 xd 3 p, (43) 

where the absolute value signifies determinant and n a are timelike vectors normal to the spacelike hypersurface. This 
is just a fancy way of writing physical momentum and space volume for a fixed-time slicing. The free-streaming 
Boltzmann equation for evolving AT is given by 

P^f - p^ c T^d p ,J = 0, (44) 

with calculated using the metric of Eq. (|31|) . After standard manipulations, one finds the temperature field 
evolution to obey 

In integrating this equation for the temperature field today, the contribution of the right-hand side of this equation is 
called the integrated Sachs- Wolfe effect, and with the right-hand side neglected, one obtains what is simply referred 
to as the Sachs- Wolfe effect. 
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If AT/T + $ is not vanishing initially, to leading approximation we can neglect the integrated Sachs- Wolfe effect 
term 2do$>, and the relevant Boltzmann equation is 



\Pc\ l ) { T 



0. 



Let us consider the initial condition 



AT 



AT 



$1 



(46) 



(47) 



for this Boltzman equation where we again remind the reader that i corresponds to the last scattering surface time 
and / corresponds to today. Inflationary computation gives us values for AT/T + $ during radiation domination 
when the modes of interest had wavelengths far outside of the horizon. We shall denote these initial conditions with a 
subscript p (primordial) such that for example $| p denotes the value during radiation domination. Starting from these 
primordial values, we would like to derive the initial condition Eq. I|47(l in terms of the "gauge invariant" curvature 
perturbation £ and the isocurvature perturbation Ss (see appendix for its definition) evaluated at the last scattering 
surface. 3 To accomplish this, we need to express everything in terms of the set of field variables {i5 7 , Sx} and then 
express the final result in terms of £ and Ss- To start off, the temperature variable is easy to exchange in terms of 
the photon overdensity since by definition 



AT 



(48) 



From now on, we will go to the Fourier space (spatially flat e lk ' x basis) and assume all of our variables now represent 
amplitudes in Fourier space. We have according to T a . a = (Eq. iJSSJl) that 



d = <L(t, k) -4$(r,/c) 



(49) 



approximately a constant on large length scales. Because the curvature invariant £ computed through the usual 



inflationary formalism is C\ = — $| 
radiation dominated era, we can write 



( / 4 where the subscript p denotes the quantities are evaluated in the 



CI 



(50) 



'p 4 

Since £ is an adiabatic invariant, £ will continue to equal C\/A even after matter domination. By Eq. I)49|l . we find 

6 



c 



_1 

4 



(51) 



for all time in the long-wavelength limit. 

Inserting Eqs. (|48|l and (|51|l into Eq. I|47|) . we arrive at 



AT 



$1 



$1 



C + 2$l 



(52) 



Note that we have here assumed that the wavelengths of interest are sufficiently long such that radiation era quantities 
are valid at the last scattering surface which is assumed to be after matter domination. 

Now, we want to reexpress $ this in terms of the isocurvature perturbations. As discussed in the appendix, the 
isocurvature perturbation of interest is 



5s = 6 



x 



—6~, 



(53) 



Hence, we need to eliminate Sx and <5 7 in terms of S s and £. 



The "gauge invariant" curvature perturbations are usually described by £ or 1Z. which are equivalent when the scales are far outside of 
the horizon. The relationship between the two arc £ = — & + 8p/3(p + P) = 1Z+ (k/aH) <t>2p/9(p + P) and we will choose to work with 

C- 
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To relate 5x with <I>, we can use the energy conservation equation (Eq. ty'33\i ) in the long- wavelength limit, and 
restricting to the non decaying mode, to write 

5 X -3® = C 2 , (54) 

where C2 is a constant. Now, since during matter domination, the 00 component of Einstein's equation (Eq. (021) in 
the long-wavelength limit is 

we can use Eq. (|54() to solve for the growing solution, which is 

(56) 

During matter domination, we thus find using Eqs. (|56|l and l|51|) that 

5s = -2 3| ( - 3 (C + $| J = -5 % - 3C (57) 

This is a remarkable feature of the isocurvature perturbations in which even if £ = 0, the isocurvature perturbations 
grow into potential perturbations during matter domination. 

We say that the isocurvature perturbations grow into potential perturbations because during radiation domination, 
if ( = 0, then the gravitational perturbations approximately vanish. Let us see how this happens in detail. Solving 
the 00 part of the Einstein equation (Eq. (|42[l ). we find during radiation domination (RD) that 

*\ P (r,k) = ^, (58) 

where as before, the subscript p denotes the radiation domination era when the wavelengths of interest are far outside 
the horizon. This and Eq. (|50() imply that in the absence of curvature perturbations £, the gravitational potential $ is 
during radiation domination. However, during matter domination, the potential <i> grows due to the existence of the 
isocurvature perturbations. One can also understand from this (the fact that during radiation domination $ = <5 7 = 
if f = 0) that if matter domination never occurs then we would have in a hypothetical radiation dominated last 
scattering surface that 



AT 
~~f~ 



, . , + hypothetical = ( 59 ) 
hypothetical 



which is in accord with the intuition that the superheavy dark matter perturbations are irrelevant if their density is 
too small. 

Combining Eq. (|57J) with Eq. (|52|l . we arrive at the desired expression 



AT 



+ *li = -^C-ffc. (60) 



C. Ci characterization of CMB 

Now consider the computation of Ci which is defined as 

(a; m a;* TO /) = Su'S mm 'Ci, (61) 

where 

AT(r f ,x,p) ^ 

^ = 2^aim{Tf,x)Yi m {p) (62) 

l,m 



and p is the direction vector of the photon. Using 

_d 

(2tt 



rj) I /',/<> Yr m {p) AT{T ^ ,P) , (63) 
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where f2 is the solid angle for p, we write the ensemble average for Ci as 



\0>lma>Vm' I 



d 3 hd 3 k 2 
(2tt) 6 



AT(r f ,k 1 ,p 1 )AT*(r f ,k 2 ,p 2 ) 



T 



T 



YiMYl'm'fa), (64) 



where we have gone to spatial Fourier space. Now, use the property that after the last scattering surface the Boltzmann 
equation Eq. (|46[) for long wavelengths is approximately 



{d Q +fd t ] 



AT 



(r 



to write 



AT(t, x,p) 
T 



+ $(t,x) 



d 3 k ( AT (it) 



T 



$ « 



$(fc) 



gik-x ^ — ip-k(r — Ti) 



Neglecting <J>(r/,:r) (today) gives 



AT(T f Xpi) „ ( AT(fc) 



T 



-i§x-k{Tf—Ti) 



Hence, we find 

(ai m a*, m ,) 



1 /" d 3 h d 3 k 2 



V J (2tt) 3 (2tt) 3 



^ AT(h) 



+ $(fc x 



AT*(k) 



T 



To leading order, the power temperature correlation function can be written as 



( AT{k x ) 



T 



$(fci 



AT*{k 2 ) 



T 



§*{k 2 



^P(k 1 )(2^f5^{k 1 -k 2 ). 



Inserting 



e - ifi a(T f -n) = £Y(21 + l^d^r/ - r,)|)fl(-fc ' k) 



1=0 



where Pi are Legender polynomials, we find 

d 3 k 



f d A k f f 9tt z 

oo 

x i' a (2fe + l)ji 2 (|fc(r/ - n)\)P l2 {p 2 ■ k)Y Vm ,{p 2 ). 



Now using the identity 



4tt 



fl HH-*) = 21 + 1 ^ 



where YJm are orthonormal spherical harmonics, we find 



{ai m a v 



dk 



6u>6 mm ,4n I ^-P{k)[M\k(T f -n)\)} 2 . 



(65) 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 



Thus, we see that Ci is essentially the power spectrum on long wavelengths up to angular projection effects of the 
Bessel function. Defining the photon travel distance from the last scattering surface 



L = T - T de 



1 



dz 



H a J Q ^ m (z + l) 3 + «A' 



(74) 
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we arrive at the desired expression for Ci as 

/dk 
—P{k)[ 3l (kL)f. (75) 

To evaluate Ci from this expression, we merely need to compute P(k), which from Eqs. I|69(l and (|60J) is seen to be 

P(k) = ^P c (k)\ l + ^P s (k)\ t , (76) 

where the subscript i denotes last scattering surface. 

To compute P$ at the last scattering surface, from Eqs. 149|) and (|54|1 we find 

Ss = C 2 - |d. (77) 

In other words, isocurvature perturbations are approximately constant far outside of the horizon, whether the quanti- 
ties are evaluated in matter or radiation domination era. Therefore, since we can express S 7 in terms of £ p in radiation 
domination using Eqs. J35J|, (|5U|l . and as 



it is convenient to compute P$ in the radiation domination era using the relationship 

Ss\ p = S x \ p - C\ p , (79) 

where one should understand the existence of nonzero Ss \ p even when Sx\ P — from the fact that the dark matter 
number density would then not trace the radiation number density. Hence, P((k) and Ps(k) are 



Ps(k)\i = 


P Sx \ p 


+ Pd P 


-Pix\ p -Pxc\ p 


(80) 


PSx = 


k 3 I 

2^ J 


d 3 r e" 




(81) 


P ( (k) = 


k 3 I 
2tt2 J 


d 3 r e" 




(82) 


Pcx(k) = 


k 3 1 


d 3 r e" 




(83) 



where all the functions with subscript p are evaluated at the radiation-dominated era. For the usual thermal CDM 
scenario, we have 

S X = h„ (84) 

which gives a zero power for the isocurvature as expected. For the isocurvature of superheavy dark matter, the cross 
correlation terms Pqx and Pxc vanish. The first term of Eq. (|76|l is the usual adiabatic contribution to the CMB 
which can be easily computed for slow roll inflationary scenarios in a standard way. The second term of Eq. I|76|l (or 
more specifically Eq. (JS^) is what we are primarily concerned with in this paper, and we will turn to its computation 
in the next section. 

Before beginning the computation of the matter correlation function, note that Eq. I|78|) allows us to rewrite the 
familiar Eq. (|76|l as 

P(k) = ±P c + ^P 6x , (85) 
where we have dropped the subscript for brevity. 



IV. ISOCURVATURE POWER SPECTRUM 



In this section, we proceed to compute Pg x (k) needed in Eq. I|85|l for the computation of Ci using Eq. (|75|) . As we 
discussed in the introduction, we restrict ourselves to a massive scalar field X minimally coupled to gravity without 
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any other interactions for X. The quantum fluctuations of the energy density during inflation is usually computed 
by using the curvature perturbation which remains constant far outside the horizon. However, the usual computation 
procedure which attempts to compute the quantum fluctuations induced about a classical field background does not 
apply to our scenario since the dark matter particles do not arise from a classical field background but from quantum 
particle production. Liddle and Mazumdar [41| had a similar cosmological scenario as the one of interest in this paper, 
but instead of real particle production, their homogeneous dark matter density was implicitly from the vacuum energy 
contribution. This difference can be readily seen by computing the vacuum expectation value of the stress energy 
tensor for a massive scalar field X with mass mx as 



3 x (T 00 (x)) = I d 3 k a 3 (A g At) 



\X k (t)\ 2 + (k 2 /a 2 +m 2 x )\X k (t)\ 



(86) 



where the vacuum state |0) is defined such that ^4^|0) = while the fields are defined by a Bogoliubov rotated boundary 
condition. The vacuum of Ref. p~H corresponds to the case with j3 k = 0, and their energy density came from the 



\m 2 \X h 



term in the integral, which is part of the zero point energy term. 



Given that the zero point energy must be taken care of by the cosmological constant problem solution, whether 
the zero point energy should be counted as an unambiguous production of particles is unclear. Indeed, in Minkowski 
space, we usually discard the zero point energy. One might argue that the zero point energy is what is contributing 
to the usual computation of generating density fluctuations in inflationary cosmology. However, this is not true. As 
we stated previously, the fluctuations that are being computed in the usual density perturbation computations are 
not the fluctuations about the classically zero energy state but about a quasi-de Sitter background of positive energy. 
Indeed, in computing the density fluctuations, the zero point energy of the inflaton is always subtracted. Hence, since 
treating the zero point energy as real particle production without specifying assumptions about the solution to the 
cosmological constant problem is speculative, we will in our computations throw away the zero point energy just as 
in Minkowski space. This is one of several important differences between Liddle and Mazumdar |41| and our paper. 

Assuming that the inflaton energy density dominates over the dark matter energy density during inflation, the 
curvature perturbations will be computed in the usual manner. This is an excellent assumption since the dark 
matter energy density to inflaton energy density during inflation needs to be much smaller than 10~ 10 to obtain 
the phenomenologically acceptable amount of dark matter today. Furthermore, because the superheavy dark matter 
field is assumed to be sitting at the minimum of its potential with zero vacuum energy, there is no vacuum energy 
contribution coming from the superheavy dark matter sector unlike the inflaton sector. In other words, the energy 
perturbation dp coming from inflation is linear in the inflaton field fluctuation S(f> while the energy perturbation 
coming from the superheavy dark matter is quadratic in X. Hence, to linear order, there is no mixing between the 
inflaton energy density and the superheavy dark matter energy density. 
The dark matter energy density fluctuation is defined as 



Sp x (x) =:T£\g):-{:T£> :) 



(A") 



(87) 



where the normal ordering is with respect to the Bogoliubov rotated operators to eliminate the vacuum energy Hence, 
we can write the correleation function as 

(6p x (x)6p x (y)) = (: T^\x) : : (y) :) - (: :) 2 . 

Explicitly, the Bogoliubov rotation is given by 

and the normal ordered stress tensor vacuum expectation value is given as 



(89) 



Px(t) 



(:T, 



(X) ,\ 



|^ {SR [a^ (xl(t)+wlXl(tj)] + |/%J 2 (\X kl (t)\ 2 + <|A% W| 2 )} , (90) 



where [dk] = d 3 k/(2n) 3 ^ 2 and we have used the fact that X^ = X_^ and defined w 2 — k 2 /a 2 + m 2 . Note that 
this vanishes in the limit 0t — > where there is no particle production consistent with our previous discussion of 
comparing our paper with [41|. Hence, the dark matter correlation function is 



(:T ( X) (x)::T ^(y):) 



(5 x (x)5 x (y)) 



(X) 



P 2 x 



- 1. 



(91) 
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With the mode mixing given by Eq. Q89[) . the stress tensor correlator can be written down straightforwardly 

x ( X kl (t)X k2 (t) - th^x kl (t)X k2 (t)+ m 2 X kl (t)X k2 (t)) e ^+^- s 



+2 f a* t a t At A t + 8 ,- ar A ,- A r + a* r 8* t At A 1 - + (3 t 8* t A t A f -. ) 

\ fei k 2 kl k 2 l^-ki k 2 -ki k 2 kr -fc 2 fei -fe 2 fea -Al -fc 2 / 



x ( (t)X k2 (t) + ^ii; (<)X fc2 (t) + m 2 A^ (f)X fa (t) 



+ f a* r a*? At At + 8 e a* r A g At + a& /3 t At A +(1 f) 1 f 4 e 4 r ) 
V fei fe 2 fei fe 2 fc i fe2 - fe i fe 2 fei ' k 2 kl -fe 2 f-fei^ -fe 2 -fei -fe 2 y 



-X* ki (t)X* k2 (t) + m 2 X* ki (t)X* k2 (t) 



-i{k\-\-k,2)-x 



[dk 3 ][dk 4 ] 



( at at At At + 8* r a t A f -. At + a r 8* r At AK + 8* r 8* r AK aK 

y fe 3 fe 4 fc 3 ki f- -k 3 ki _ fc3 ki fc 3 ^-fe 4 fe'3 - ki ^-fe 3 ' -ki -k s -fe, 

x (x k3 (t)X ki (t) - ^^X k3 (t)Xki(t) + m 2 X k3 (t)Xki(t)^ e'^+^-v 
+2 f at at At A % +(3 a % A E At + a} (3* t At A f r + 8 s (3* r A t ^ r ) 

V ^3 ki k 3 ki ^-k 3 ki -fe 3 ki k^' -ki fe 3 -fc 4 ^—k 3 r_ ki -k 3 - ki ) 

x ( X* k3 (t)X k4 (t) + k -^Xl (t)X ki (t) + m 2 X* k3 (t)X ki (t)) e-^-^ 



+ (a* r a* r At At +8 t a* r A t At + a*, [3 t At A ,- + 8 ,- [3 t A ,- A A 

\ k 3 ki k 3 k 4 ' -k 3 ^ -fc 3 ki k^-ki k 3 -ki ' -k^-ki -fc 3 -fe 4 J 



x ( X* k3 (t)X* k4 (t) - ^X* k3 (t)X* ki (t) + m 2 X* k3 (t)X* ki (t) ) e -*(S»+Sd-» 



(92) 



Although not rigorously justified, we can approximate to within an order of magnitude \8 k \ <C \a k \- In the expansion 
with f3 k — > 0, the correlator is simply 



:T o W(z)::T o W(0 



The power spectrum is 
fc 3 



1 /■ d z kid z k 2 



(2tt) 6 



|afei| 2 |afe 2 



^ii^fcj — y~ ~X kl X k2 + m 2 X kl Xk 2 



i(k 1 +k 2 )-(x-y) 



+0(8 k ). 



(93) 



P 'x(k) = ^ I d 3 re~^(5 x (x)6 x (y)) = 



1 fc 3 1 



+a- 2 (l^l 2 - costf) X kl X V k2+k2 _ 2kikcosf) 

Note that the negative unity in Eq. (|91|l disappears because 

A: 3 



dk\k\ I dcos9 



x k ,x 



m 2 X k ,X 



1 v /fe 2 +fcf-2feife( 
2 



^/fc 2 +fef-2feife COS 



2tt 2 



0. 



In the flat space limit, we find 



• ^00 ( x ) ■ • -^oo \ x ) 



If we neglect the kinetic terms, we find 



1 f d 3 k 1 d 3 k 2 
2 



(2tt)6 (2E kl )(2E k2 ) 



-Eh, E, 



ki ■ k 2 



kx^k 2 



m 



1 TTL 

(2,)a pTO 'TWW)' 



(94) 



(95) 



(96) 



(97) 
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which is what we expect since in that limit p ~ m 2 X 2 and 

(: X\x) : : X\y) :> = 2{X{x)X{y)) 2 . 

Note that in general, the kinetic part of the correlation function cannot be neglected. 

To evaluate the correlation function in the curved spacetime of interest, we must compute the wave function, which 
solves the following Klein Gordon differential equation: 

X'{(t) + 2HX' k (r) + (k 2 + a 2 (T)m) X k (r) = (99) 

Note that, since the Dark Matter field has no VEV, the equation of motion does not couple to the metric fluctuations 
at first order. 

We write the evolution of the scale factor in conformal time r as 

-^H = -eH 2 , (100) 
a dr 

where H = a~ 2 da/dT, and e is the usual slow roll parameter (potential derivatives) which we take to be approximately 
constant (This is true at first order in the slow roll parameters). With the boundary condition that a(r,) = a% — 
— l/(HjTi), we find 



1 



Hit 



1 + e 



1-Ii -Inl 



(101) 



where one recognizes the usual dS scale factor —1/(Hjt) in the limit that e — > 0. However, using this would make 
an analytic solution impossible since a 2 will have terms as well as ln(r/ri) terms while a" /a will have 1/r 3 terms. 
Hence, we will use for the dark matter correlation function in the dS approximation. 

Still, we can approximately take into account the changing H for the isocurvature perturbations during slow roll 
inflation as follows. Since the amplitude of the perturbations are approximately frozen when the physical wavelength 
crosses the horizon, we will assume that H relevant for the density perturbation to be k dependent such that 

(102) 



a{t k )H k 

where t k corresponds to the time at which k/a = H k - Now, solving the slow roll equation 

H 



JJ2 '■ ( 103 ) 



we find 

1 1 



„ -<tk-ti), (104) 



where Hi corresponds to the expansion rate at some initial time ti. Using approximate dS expansion 

1 , (a k 



(t k - tl ) = ±-hx(^\, (105) 



we write 



or equivalently 



-J- + -J- = -^f-l. ( 106 ) 



Hk Hi Hi \ai 



• • (^) (1H7) 



Using Eq. (|102(l . we arrive at 



H k = Hi[ — — ) + 0(e 2 ). (108) 
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Hence, the mode function can be approximated as 



Xh — 



2cPl 2 JTT k 



v \aH k J> 



where 



V 4 \H k 

The mode functions satisfy the usual normalization 
In the long- wavelength limit, we find 



.(-_ 1 )3/4 2 -l-^ e -W7r/2 e -H fc t(3/2- l /) fl -i'-l/2 fe - 



-2vH k t 



HI 



where we have used a = e k . When v is real and positive, this behaves as 



X k « _4=(-l) 3 /*2- 1 - , 'e <w / a f -4 



4' 7 



(109) 



(110) 



(111) 



(112) 



(113) 



which can become large when k/(aH k ) <C 1 and v ^ 0. When ^ is imaginary and positive, the wave function behaves 
as 



-1 



X k » -= (-l) J / 4 2 



7T/2 / fc 



2 | V | 



1 



V aH k J \/H k w 



■T(-i\v\), 



(114) 



whose magnitude is essentially independent of k, as it can be understood by noting that this limit corresponds to 
to 3> k. 

Let's see how the mode function X k scales after inflation ends. Assuming that the scale factor scales as 



we find for the mode equation 



/- 



f = o, 



(115) 



(116) 



where X k (t) — a 3 / 2 /(i). Neglecting the k term in the long- wavelength limit and defining £ g = 3/2q — 9/4, we find 
that in the limit to 2 3> H 2 , 



X k 



,3/2 



while in the limit m 2 <C H 2 that 



Xu 



7,3/2 



i (l + x /l-4? 2 « <I )/2_ 



(117) 



(118) 



Hence, we see that since H ~ eventually, to 2 3> -ff 2 and X k will fall like a 3 / 2 . Hence, for any quantity we 
compute involving X k mode function in dS space, as far as the scaling with a is concerned in the absence of further 
interactions, we should freeze its value at the point when to = H and then scale it as a~ 3 / 2 . 

Because of the k suppression, the term a~ 2 (|fci| 2 — |fci||fc| cos9)X kl X ^/ k 2 +k 2_ 2kikcoeg in Eq. lf9"Hl should be negligible 

unless there is a cancellation of the other terms in the limit that |fc| — > 0. Quantitatively, the ratio 



1 y / k 2 +kj~2k 1 kcost 



lim 

fc,fci^o m z X kl X 



y / k 2 +kj-2k 1 kcost 



m 2 

4to 2 



(119) 
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Hence, as long as m 2 ^ 9i? 2 /4, there should be no cancellation. Hence, the power spectrum Eq. I|94(l becomes 



1 fc 3 1 
— o~— ,„ .„ / dkiki I dcos( 
p\ 2tt (2tt) 3 ' 1 ' 



Xu. X 



A //£ 2 +/cf-2/ci/ccos0 



m 2 X fel X 



1 fc 3 1 



dk%ki / dcos( 



^/fc 2 +fcj-2fcifccos0 



1 \A 2 +fcf-2feifc 



^/fc 2 +fc 2 -2fcifecose 



-2m 2 Re 



Lfc l^^fc2 +fc 2_ 2fclA . cose ^l^ fe 2 +fc 2_ 2felfecos6 , 

In the language of j^J, we have 

P x (fc) = ^l|X fc | 2 , 

27T Z 

which implies, keeping only the term proportional to m 4 , that 



(120) 



(121) 



P Sx (k) = 



k 3 1 



Px 2tt (2tt) 



dfcxfc 2 / dxP x {kx)Px{\ki-k\) 



2tt 2 2tt 2 



= ™k 3 — I d 3 h 



k\ \ki~k\ 3 P 2 x 8tt 



Px(ki)P x (\h-k\) 



k 3 \h-k\ 3 



(122) 

which matches the result of Ref. 41]. However, unlike their analysis, we do not neglect terms involving 
Xh-, X / ,,,,, — zrr—, 7, which can be a priori of same order of magnitude as the m 2 Xk, X / ,„,,, — jrr - ; 7, term they 

1 y k 2 + k( — 2fcife cos 8 ^ o in y'k 2 +kf—2k 1 kcos8 J 

accounted for. We will see the resulting correction is close to a factor of 2. 

We can express the isocurvature power spectrum Eq. I|12Q(I in terms of the power spectrum of the dark matter field 
as follows: 



P Sx (k) = 



k 3 1 



d 3 k x - 



Px^J k 3 \k-h\ 



Px(h)\ +^/ a 6(2vr 2 ) z J \\Px(\h-k\)\ +\k-k 1 \ya e (2ir 
WPxik^Pxdk-hl) 



+ m 4 P x (k 1 )P x (\k 1 - k\) + ^ (Px(fci)Px(|fc - fell) - k 3 \k- k 1 \ 3 /a 6 (2n 2 ) 2 ) L . (123) 



Note the presence of additional terms with respect to Eq. (|122|l due to the kinetic contributions to the stress tensor 
appear as derivatives of the power spectrum. 

At the end of inflation, the Hubble expansion rate is much smaller than the value at which the isocurvature 
perturbations are generated. For example, in the chaotic inflationary scenario that we consider, we have 



(124) 



where Hi is the expansion rate during the time that the isocurvature perturbations are generated. Since gravitational 
particle production is determined by 



mx mx V50 

in ~ ^ ' 



(125) 



and since to have enough particle production for dark matter with Trh P$ 10 9 GeV, 4 we must satisfy the condition 

m x 



H, 



<7, 



we must have 



m x 
Hi 



< 1. 



(126) 



(127) 



4 We will see later that the isocurvature bound of interest is within the regime where Tjuj < 10 9 GeV. 
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This implies that v is real and positive, 5 which means that the wave function will behave as 

x k » _- Lf-i) 3 /^- 1 -"^/ 2 ( -4 ^ J " 

Vtt i ..II. 

and the spectrum 



\aHkJ ^/H^ 3 



r(i/), 



7,3 



-2v 



where 



o2i/-3 

a^— Hi»| 2 . 



(128) 



(129) 



(130) 



When we take the time derivative of Px(k), we will only account for the time dependence in a since that is the main 
time dependence information contained in the approximate Bessel function equation. Hence, we find 



P x « A 



k 3 ( k 



Hka 3 \aHk 



-2v 



(2u - 3)H k 



(131) 



where in the second equality, we have assumed approximate dS expansion again. 
Finally, using Eq. (|108|l . we find 



Px(k) = A 



k 3 



PX = A^r 



if/a 3 \aH! 



k 3 ( k 



a 3 \aHi 



(2u-3) 



where 

n = e - 2i/(l + e). 

Hence, we can write the simplified form of the power spectrum as 



(132) 



(133) 



P5x{k) = ! k\L dUu2 \ 256aW Hl {-Jk 



-(l+2n) ( / U \ 2 (™- 

(u)-M 1+4A 5 ' 



\aH! 



i 2(n- £)(3 _ 2l/) 2 7r 4 



1-u 



-n+2 



+ \l + u\ 2 ~ n 
+A 



-4A»(3-2i,)V f ft \ 2( ^ £) |1 _ u|2 („V 

| 1 + u |2(n-e) 



i 2(l-e) + n \aHj 

1 4A 2 (3-2z/) 2 tt 4 / fc ^ 2(n ~ £) 



\2-n 2(l-e)+n \ a ^/ 

2n— 1 4. 

[|i + u |"+ 2 -|i- u |"+ 2 ] 



i 7 H 3 (n + 2) 



m 
~4a? 



-1 -2A 2 (3-2i/) 2 / fc 
+ u(2- e + n) \^Hj 



i + «r e 1 -u\~ e u r ' 



x (|1 - U |"+ 2 |l + u| £ - |1 - U | £ |l + u\ n+2 ) 



(134) 



One can easily check that the term proportional to m 4 matches the integral expression of |4l| • 



5 With nonminimal coupling, v may easily be imaginary in which case Eq. 11141 will be relevant instead of Eq. 11131 . In that case, the 
spectral index for the isocurvature perturbations will not be as sensitive to the mass mx since ix)\ exp(x7r/2) scales as X/s/x for 
real x > 1. 
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TABLE I: Coefficients to the integrals governing the order of magnitude of the power spectrum is presented. The estimate in 
the third column includes the contribution from the integral multiplying the coefficients in the expansion of the power spectrum. 
As explained in the text, the divergent integrals were cut off at fc max = a e H e and the k value of the spectrum set at k — aiHiy 
where <ij is the scale factor at the "beginning" (near 50 efolds from the end) of the inflation and a e is the scale factor at the end 
of inflation. The degree of divergence in the integrals is given in the fourth column. It is clear that (ah contributions to the 
power spectrum dominate because these terms are not diluted by the enormous scale factor. For C2I2 (where I2 is the integral), 
the divergence is one degree smaller than the degree of divergence in each term in the integrand because of cancellation. 







order of magnitude 


degree of 


coeff. 


value 


contribution including integral 


divergence 


Cl 


(256A 2 Tr s )' 1 (k/aH I ) 12 - 2n ' 


6 x 10- i (a e /a) 10 (a i /a) 2 (H e /H I ) w y 2 


10 


C2 


4(3 - 2v) 2 n 4 (ni - 1 - e)- 1 (256w s y 1 (k/aHi) 6 - 2e 


7 x lQ- 5 {a l /a) 3 {a e /a) 3 {H £ /H 1 f{m x /H I fy 3 


3 


C3 


ci/(5 - ni) 


6 x lQ- 4 (ae/a) 10 (ai/a) 2 (H e /Hi) w y 2 


10 


C4 


(rat — 1 - e)c 2 


-7 x 10- 5 (a 8 /a) 4 (a e /a) 2 (H e / J f/ / ) 2 (mx/H/)V 


2 


C5 


-(3 - 2^) 2 (647r 4 )- 1 (n i - 5))- 1 (k/aH I f- 2c 


10- 5 (a l /a) 4 (a e /a) 2 (// e / J ff / ) 2 (mx/// / )V 


2 


C6 


-m 2 x {A-K 4 'H'j)- 1 {k/aHi) & 


-3 x W-^ae/afiai/afiHe/Hifimx/Hjfy 3 


3 


ai 


A 2 {?> - 2^) 4 16- 1 (1 + 2e - my^k/aHj) 2 ^- 4 ' 


8 x l(r 6 (mx/#/) 8 




02 


A 2 (m x /Hj) 4 (l - mr'ik/aHj) 2 "' 


6 x 10- i (m x /H I ) i 




O3 


A 2 (3 - 2u) 2 {m x /H I ) 2 2-\l - m + e)- x {k/aHi) 2 ^-^ 


l(T 4 (mx/H/) 6 





TABLE II: Integrals multiplying the coefficients in Table 



integrand for the measure f Q max du 

T 4— n; I I I 5 — Tlr A 

I 2 tZ-^fll + ul"'- 1 - 26 - |l -u\ n i- x - 2c ) 

h u 4 - n '\l + uf- n ' 

h u 2(ni - £ - 2) |l - u| 5 - ni 4 

h u 2(ni - 2 - £) |l + u| 5 - ni 4 

h u 2 

F 1 u n '- 2 - 2 '(\l - u\ n '- 1 ' 2 ' L - 1 +u\ ni ~ 1 - 2t ) 

F 2 u n '-' 2 (\l-u\ ni - 1 -ll + u]^- 1 ) 

F 3 u n '- 2 -'\l + u\- e \l-u\- e (\l-u\ n '- 1 \l + u\ a - |1 - ^ni + ^l"'- 1 ) 



Since we know from Eq. I|133|) that the powerlaw index n is generically a negative number close to —3 (since m/ Hi 
is expected to be much less than 1), we will define a new power law index 



The u integrals can be separated as 



ni = n + 3. (135) 



(136) 




where all the coefficients are shown in Tableland the integrands to the integrals are shown in Table HT1 The estimate 
in the third column is given with any divergent integrals cut off at fc ma x = a e H e corresponding to it ma x = a e H e /k and 
k = aiHiy where y < 10 4 is a scaling parameter relevant for the CMB: e.g., 

du [u 4 ~ n > (\1 + ul"'- 1 - 2 '- ~ \1 - ul 71 '- 1 - 2 ')] . (137) 

This is a sensible cutoff since the power law form of the wave function and the decoherence of quantum fluctuations 
no longer holds when k > aH , and the shortest wavelength that gets stretched beyond the horizon during inflation 
is k — a e H e . It is clear that aili contributions to the power spectrum dominate because these terms are not diluted 
by the enormous scale factor (a e /aj > 10 22 ). This means that the prediction is insensitive to the cutoff and the 
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FIG. 1: Plots of the functions important for the isocurvature spectrum. Here, e = 0.01 for F\ and F3. 



prediction is robust. Note that in making the order of magnitude estimate in the the third column, we have assumed 
that m ~ 0(m 2 x /H 2 ). 

The spectrum can thus be written simply as 



H 8 1 



A 2 (Z-2v) i 



2n t -4e 



16(1 + 2e-m) \aHi 
A 2 (3-2v) 2 (m x /H I ) 2 ( k 



2(1 -m+e) 



aHi 



F 1 (n l ) + 

2(nj-e) 

F 3 (ni) 



A 2 (m x /H I ) 4 ( k 



ni 



\aH! 



2 n 1 



F 2 {ni) 



(138) 



Strictly speaking, this expression unfortunately cannot be used beyond the end of inflation since the wave function we 
have used to compute the spectrum is no longer valid even approximately. Indeed, as we discussed earlier, when the 
Hubble expansion rate falls below the mass of the dark matter, the mode function scales like l/a 3//2 . However, in the 
spirit of the usual scalar density perturbations during inflation, we will simply treat the isocurvature perturbations 
as a classical fluid at the end of inflation. Owing to Eqs. H77|l and this means that for long wavelengths P$ x is 
approximately frozen at the end of inflation assuming that the final particle production density is simply scaled back 
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to the end of inflation. Hence, our final expression for the spectrum is 

Hf A 2 ( k \ 2ni ( (3-2z.) 4 / k \~ 4e , . (mx/HjY 



p|-(te) 4 Va e i?/ / Il6(l + 2e-n ; ) V a <A/ 1 - "J 

(3-2^) 2 (m x /iJ / ) 2 / k 



2(l-n/+e) V a ei?/ 
The three function Fi(n{) are shown in Figure ^ 



F 3 (n,) (139) 



V. CMB SPECTRUM 



In this section we compute numerical values for Ci contribution coming from the superheavy dark matter isocur- 

3' 



vature perturbations within the context of a specific slow roll inflationary scenario, namely that of V((f>) — ^m^cj) 2 
scenario. In particular, using Eq. (|85|l . we define the isocurvature contribution to Eq. (|75|l as 



For fi\ = 0, fi m = 1, and Zd ec = oo, £ defined by Eq. (|74|l becomes the well known expression 2/ (ao-f/o)- Numerically, 
we will instead take more realistic values of Oa = 0.7, Q m = 0.24, and Zrf ec = 1100, which yields 



Now, we will use the approximation 



where v is defined as [c.f., Eq. I|110|) ] 



a Q H ^ ^ 



m sa 3 + e-2(l + e)i/, (142) 



/9 _ /m\ 2 



4 Vff// 



(143) 



and the slow-roll parameter is evaluated at when the longest wavelengths of interest is leaving the horizon and Hi 
corresponds to the Hubble expansion rate at that same time. Using the well known formula 



k \aeHj lJ y n {a e H I L)iT(l + 2-q/2)T 2 (3/2-q/2 



we find 



_ 16 Hf 12 4 "- 6 [ (3-2^) 4 2 2 "'- 4e - 3 7r r(Z + n, - 2e)r(2 - 2n, + 4e) 



16(l + 2e-m) v J (a e H I L) 2n '- 4,i T(l + 2-ni + 2e)T 2 (3/2-ni + 2e) 



(mx/H!)*. 2 2n '~ 3 n T(l + n,)r(2 - 2m) 
-F 2 (ni)- 



1 - nj v w K#/£) 2 "< r(Z + 2 - n ; )r 2 (3/2 - m) 
(3-2v) 2 (m x /H I ) 2 2 2n '- 2e ~ 3 n T(l + m - e)Y{2 - 2n/ + 2e) 



2(1 - ni + e) ° v w (a e iJ/L) 2 "'- 2e r(Z + 2 - n t + e)r 2 (3/2 - m + e) 



(145) 



Note that Eq. (|144fl says that the effect of the Bessel function integral is to merely set k to 1/L and multiply the 
{k/{a e Hi)) q \} := i/i J by a dimensionless number of order 0.1. 



As discussed before, note that we are also effectively neglecting Bogoliubov mixing effects for the computation of the two point function. 
This is also not rigorously justified in this paper but is consistent with the usual assumptions made in computing scalar perturbations 
in inflation. 
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FIG. 2: The region of parameter space excluded by the isocurvature constraint. The shaded area bounded by dot-dashed curve 
is ruled out due to overproduction of isocurvature perturbations. The bound on C; we took was 5 x 10~ 10 corresponding 
to Cil(l + l)/(27r)|i=2 = 3532 (fiK) 2 . The dashed curve indicates the uncertainty in the isocurvature computation due to the 
uncertainty in the particle production computation. Below the solid curve, there is not enough gravitational production of 
particles to have sufficient dark matter (Qxh 2 < 0.1). Above the dotted curve, there is overproduction (Qxh 2 > 0.2) of dark 
matter. Although this specific example is for V = ^m^rf) 2 inflationary potential, as we argued in Section 2, the constraints are 
general in that superheavy dark matter scenario requires a reheating temperature of above 10 s GeV. 



Finally, for the numerical analysis, we need an expression for p e . As we explained in the appendix, the infrared 
contribution to this density is ambiguous. Although an accurate computation of the non-infrared contribution can only 
be computed on a model by model basis, the mass bound coming from the isocurvature contribution is not sensitive to 
the uncertainty because of the large number in the exponential. On the other hand, the uncertainty on the reheating 
temperature is the same as the one in the density of particles produced. Hence, the reheating temperature bound is 
uncertain by a factor of about 10. 

For the slow-roll inflationary scenario of = ^m 2 ^ 2 , one can compute the non-infrared contributions to the px(t e ) 
to be bounded as shown in Eq. (|14[). For a more general class of slow roll inflationary models, the particle energy 
density can be estimated by Eq. I|15fl . To obtain the numerical value of a e HjL appearing in Eq. (|145fl . we use the 
standard reheating relationships 

/ \ 1/3 /rp \ 4/3 

a e PRH\ ,„ „ „ 1n -62\l/3 / J -RH 



(2.2 x KT 62 ) 1 / 3 ('- 

1 ' V 



aRH V Pe J \GeV 

a RH ( g.8(h) \ 1/3 T _ 81xlQ _ 14 GeV (W6) 

CDM 

in the V = \m 2 ^4> 2 slow roll inflationary model. The somewhat arbitrary bound on we took was 5 x 10~ 10 

corresponding to C;i(Z + l)/(27r)|/ =2 = 3532 (pK) 2 in units of temperature. This rather large upper bound was taken 
to obtain a conservative bound on the mass of the superheavy dark matter. Because of the largeness of a e HjL in Eq. 
(JHSJl, the mass bound is not sensitive to the exact bound of C\ we use. As one can see, since the dark matter must 



a o \9*S\tRH) J Tub Trh 
In Fig. |21 we plot the parameter space allowed by the combined constraints of isocurvature fluctuations and Q 
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compose the CDM of the universe, we have a robust bound on the reheating temperature of 

T RH > 10 8 GeV. (147) 
and on the mass of the superheavy dark matter of about 

it- 6 ' (M8) 

As argued in Sec. 2, this bound on the mass and the reheating temperature should be robust for most slow roll 
inflationary models. The rheating temperature bound has an uncertainty of about a factor of 10 and the mx / H e 
bound has an uncertainty of less than unity. Note that one cannot naively extrapolate using our results to situations 
to large mx/H e (corresponding to small lower bound on reheating temperature) because there, the gravitational 
production is insufficient to compose the dark matter and the computations have assumed that the gravitationally 
produced superheavy dark matter dominates the matter density during matter domination. When the gravitationally 
produced contribution becomes negligible fraction of CDM, we expect the isocurvature bounds to disappear. 7 Unlike 
in section 2, the tensor to scalar amplitude limit does not play much role here because the tensor to scalar ratio here 
is predicted in this specific inflationary scenario to be much less than unity (of order lOe ss 0.1) on long wavelengths, 
consistent with r$ = 0.1. It is important to note that all of the results in this section are in agreement with the more 
general discussion of Sec. 2. 

One might worry that since mx/H e considered in Fig. [21 corresponds to a number larger than 1, we may have 
violated our assumption that v is real. In most inflationary scenarios, this constraint is not violated since Hi > ZH e 
even when slowly rolling. (For example, Hj > 7H e for V = \m?^ 2 type potential.) 

Note also that numerically, the relative contributions to C\ X ' of the F-z, F3} terms are 

(term Fi) : (term F 2 ) : (term F 3 ) = 1 : 11 : 7 (149) 

at mx/H e — 6 and Trh = 10 9 . These ratios are typical in the region of interest. Hence, the kinetic contribution 
gives almost a factor of 2 correction to the isocurvature perturbation computation. 



VI. CONCLUSION 



Isocurvature perturbations are generic prediction of inflationary cosmology. We have explored the isocurvature 
perturbation constraint on models of superheavy dark matter minimally coupled to gravity. Surprisingly, there is a 
robust lower bound on the reheating temperature for these models of about 10 7 (0.2/rs) GeV where r$ is the bound 
on tensor to scalar. This means that, if the superheavy dark matter scenario is embedded in supergravity models 
with gravity mediated SUSY breaking, the gravitino bound, which is typically an upper bound on the reheating 
temperature of about 10 8 GeV, strongly squeezes these models. If, for example, LHC data favors a gravity mediated 
SUSY breaking scenario in which the cascade of gravitino decay reactions favors an upper bound on the reheating 
temperature of 10 7 GeV, then superheavy dark matter scenario can be ruled out. Note also that because of the 
presence of tensor to scalar power limit rs, improved experimental bounds on this quantity can also squeeze this dark 
matter scenario. 

There is also a corresponding mx / H e > 0(5) bound on the mass mx of the superheavy dark matter. This 
implies that the gravitational particle production is in the exponentially suppressed regime since the exponent of the 
suppression is of order —2mx/F[ e : i.e. noninfrared modes make up most of the dark matter. Moreover, as can be seen 
in Eq. I|29[) . one cannot lower H e to less than 3 x 10 12 GeV while keeping mx/H e fixed and Trh < 10 9 GeV to obtain 

sufficient abundance of dark matter because the overall dark matter density scales as m^/ 2 . 8 Hence, we see that for 
the superheavy CDM scenario of gravitational particle production, CMB data forcing the isocurvature contribution 
to be small favors a relatively large gravitational wave amplitude. 
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APPENDIX A: MEANING OF ISOCURVATURE PERTURBATIONS 

Given TV particle species with their energy densities different from the average energy density, the total energy 
densities contribute to the curvature perturbations (because local energy density determines the gravitational poten- 
tial). The N — 1 remaining energy density field degrees of freedom can contribute to what is often called isocurvature 
perturbations since N — 1 degrees of freedom can be assigned any value while keeping the the total energy density 
and hence the curvature constant. 

To write down explicitly the usual condition for adiabatic perturbations, denote the energy density deviation from 
the average density of particle species i as 

Si = ^- (Al) 

For perturbations to be adiabatic, we have 

5{s t a 3 ) = (A2) 

where is the entropy density and a is the scale factor. Since the photon entropy is conserved assuming adiabatic 
evolution of the photons, Eq. (|A2|) implies 



5{ Si /s 7 ) = 0. (A3) 
Since according to first law of thermodynamics, we have s 7 = (p 7 + P 7 )/T = |p 7 /T, we find 

5(sj/s 7 ) _ Ss t Ss 7 _ 5si ST _ 5si 3 _ 

7 — / — V" — — — — t ~ — A 1 ~ ^ ' 

ySijSry) Si 5 7 Si 1 Si Q 

as the condition for adiabatic perturbations. 

Now, we can derive a useful relationship in terms of equation of state. First, the first law of thermodynamics gives 

d(pia 3 ) = T^Sia 3 ) - Pida 3 + ^d(nja 3 ). (A5) 

Setting the entropy to zero gives 

dpid 3 = -d(a 3 )[pi + Pi - pint] + pidn % a 3 (A6) 



or equivalently 



S t = 3^ 
a 



Pirn _ f 1 Pi_ 
Pi 



Pi ■ (A7) 

Pi 



The adiabatic condition Eq. (|A2|) gives 



Ssi 8a 5 l -p i 5n i /pi 

Si a 1 + Pi I pi - piiiil pi 



where we have used Eq. (|A7|) . Hence, using Eqs. (|A4|) and 1A8L the adiabatic condition can be written as 

Si - pMi/pt 3 _ (A9) 



1 + Wi - Pinjpi 4 

Note that here all evolution was adiabatic. For the nonrelativistic particles, this expression simplifies further since 

Si- fiiSrii/pi 3 Siil-Hi/mi) 3 3 

- 7^7 = — I —J- 7 S 7 = °i - 1 S 1 = °- ( Al °) 



1 + Wi-piTii/pi 4 l-pi/m 
which is the usual condition that we see in the literature 
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It is interesting to rederive Eq. (|A9|I from a kinetic theory point of view. From a microcanonical ensemble con- 
struction, we can write the entropy density as 

'" " '/-•: 1 In/,:. (All) 



(2tt) 3 

where p is the comoving momentum and / is the phase space distribution function. This directly gives 

dp __ 3 , £ n Sa 



Ssi = - I j^a- i 5f l \nf l --Z—s l . (A12) 



On the other hand, the adiabaticity condition Eq. <|A2j) implies 



6si = -3— (A13) 



Hence, we are lead to conclude 

d 3 p 



(2tt): 



a- 3 6filnfi = 0. (A14) 



In the near equilibrium case without degeneracy, ln/^ < 0. If in addition, if Sfi = 0, we can conclude brii 
—3a~ 1 5a J d 3 p(2ir)~ 3 a~ 3 fi = —3a~ 1 rii5a. However, generically, the constratin is only given by Eq. l|A"H|) . 
If we assume equilibrium situation without degeneracy, we have 



d3p ;a- 3 f,(l + ^)=n l + ?l-^. (A15) 



The constraint Eq. (|A14(I gives 



Since 



we have 



2tt) 3 " " V~ ' T, 



d3 P .- 3 „ (*k 



(2tt) 3 j 'V Ti 



a~ 6 Sfi ^tf^ = 0- ( A1 6) 



a- 3 E l f l , (A17) 



a ri "7 (2tt)3 a J (2tt)3 E % jl 

-3-( Pi + + Mi / 7^S-a- 3 ^ (A18) 
a 7 (27r) ,:> a 



where we have used the adiabaticity constraint. Also, we have 

Sm= I ' ^a-Hh-^m. (Ai-n 



Hence, we write 



arriving at 



&CL ( &CL \ Set 

Spi = -3 — (pi + Pi) + [ii ( fcj + 3 — m ) = -3 — (pi + Pi - fiim) + ^iSrii, (A20) 



8 Si _ _ ^ Pi - fM^JH _ Sj - piSrij/pi (A21) 

Si a pi + Pi-pim 1 + Wi - piUi/pi' 



25 



This is precisely the formula that we derived with thermodynamics. 

Now, we can easily show that adiabatic perturbations usually imply nonvanishing curvature perturbations. Accord- 
ing to Eq. I|A7(I . we have 



8_p 
P 



- P iSi = 



Sal ^ 
-3 > pi 

a p ' 



(1 



mm 



- p-tSrii = -3—- y^TtSt + - pt8n l 
p . a p / — J p 



(A22) 



which is nonvanishing as long as piSrii term does not cancel the entropy term. If the adiabatic perturbations conserve 
particle number (which would be true for adiabatic perturbations of nonrelativistic particles such as baryons which 
carry nonncgligible chemical potential), we can write 



Sn, 



6a 
-3 — rii 
a 



(A23) 



which would imply there being no cancellation and Sp/p =^ 0. On the other hand, it is clear that one can have 
energy density perturbed without conserving entropy locally. Hence, curvature perturbations do not imply adiabatic 
perturbations. 



APPENDIX B: PARTICLE PRODUCTION IN THE LARGE MASS LIMIT 



1. General consideration of the exponential suppression 

In this appendix, using |l5j . we give an approximation of the exponential damping of particle production in the 
regime mx/H — ► oo. 

Consider the approximate |/3| 2 formula in |15| 



1/9*1" w exp ^ 

C(rj) = a 2 ( V ) 
where r is the real- value solution to the equations 



4(fc/yC(0) 2 



Amx 



6 J m x 



m xy /C"(r)/(2C*(r)) ^C" (r)/(2C 2 (r)) _ 
„ 1\ 6a" /a 3 " 



t-C"'(r)+C'(r) 
6 

^C"{r)+ k2+m x C ( r ) 



m 



(Bl) 
(B2) 

(B3) 
(B4) 



A' 



with p a pure imaginary number (that also needs to be solved along with r) . The primes are conformal time derivatives 
as usual. Note that this approximation is valid only if the conditions 



C""\r) (k 2 +m 2 x C(r)) 
C"'(r) 10m 2 x C"(r) 

\C""(r)\(k 2 + m x C(r)) 
[C"{r)Y 



6m x 



< 1 



< 1 



pf < 



(B5) 

(B6) 
(B7) 



are valid. 

To use this formula, we need to solve for r given an inflationary model with potential V((j>) which determines C . 
Solving Eqs. ljB3(l and l|B4p for r is equivalent to solving 



C'{r)C"{r) 



1 



C{r)C'"{r) ~ ' C{r)m x 
from which we see that r is a k dependent quantity in general. Now, consider the large mass expansion: 

C(r,)*a 2 ( V ) [l + 0(H 2 /m x )} . 



(B8) 



(B9) 



2G 



We then write Eq. (|B8|) as 



2a'{a' 2 + aa") 
a{ia'a" + aa 1 ") 



k 



a 2 (r)m 



(BIO) 



x 



where the primes are with respect to conformal time. Note that this equation gives r as a function of k. Hence we 
will implicitly mean r — r^. 

Using standard equation of motion relationships one can easily derive 9 



Q = 



2H[-(<j> 2 -V)/3 + 2H 2 } 
#(-2</> 2 /3 + 8U/3 + AH 2 ) + Vjj)' 



(Bll) 



Now, substituting 4> 2 = 2\3H 2 - V], wc find 



<l 



which is a remarkably simple formula. 

Substituting for (f> in terms of H and V, we find 



d 



2VHdt 



VW)) 



(B12) 



2HV 



AHV + sign(<j>)V 4> y/2[3H 2 -V] 



(B13) 



Except for the Taylor expansion approximation and the saddle point approximation that led to Eq. IjBlJI . we have 
only used one additional assumption of Eq. l|B9fl thus far. Now, we reparamctcrizc r = by writing 



^ 2 (r) « XV(^(r)), 



(B14) 



where A is a numerical coefficient of 0(1) [anticipating the fact that r will be near the end of inflation]. Again, since 
r = r*fc, we have A = Afc. The utility of this parameterization is to bound A thereby bounding the final Bogoliubov 
coefficient. 

We can then write 



1 (\ 



3 V 2 



1 V. 



(B15) 



Using sign(0) = — sign(U0) with H > 0, we find 



1 



3A \V a 



\ 2]j l + A/2 V 
Hence, we arrive at the condition determining r [Eq. (|B10|) ]as 

\vA" 



(B16) 



1 



3A 



l + A/2 V 



1 



k' 2 



a 2 (r)m 



(B17) 



x 



where the k dependence cannot be trusted for k/a < H. However, this does not matter here since we are concerned 
with H <C mx ■ 

Now, let us evaluate 



f~i/i 

W ^2CP= 2H2 



(B18) 



9 From now on, we will be using the convention M 2 ; /87r = 1. 
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which appears in the denominator of the exponential in Eq. <|B1|) . Since Einstein's equation with inflaton domination 
gives 

d=~(^-V), (B19) 

we find W — V . Hence, we find the approximate Bogoliubov coefficient amplitude to be given by Eq. IjBlfl with 
C"/2C 2 = V and r given by the solution to Eq. I|B1 7fl . Explicitly, the Bogoliubov coefficient is approximated by 

^V^ 2 - Am * ) . (B20) 



\Pk\ 2 



exp 



2. Application 

Consider the quadratic monomial inflaton potential: 

1 



V = -mi 



We find 



|/3 fe | 2 re exp • 



4^2 



(fc/qp)) 2 m£ 



where we have restored the reduced Planck mass. We use Eq. (IB 17(1 and find 

<Kr) ^ 



3A k 2 /a 2 (r)+m 2 x 



Mr 



Hence, we arrive at the Bogoliubov coefficient 



1 + X/2 2k 2 /a 2 (r) +m 2 x ' 



1 & 



exp 



2 + X k {2k 2 /a 2 (r k )+m x 



3A* 



m^mx 



(B21) 



(B22) 



(B23) 



(B24) 



Note the nontrivial factor of 2 in front of the k 2 /a 2 . 

Now, we have to solve for X k in Eq. (|B14|) . Let us look for the field value at which 

2 (r fc )~ \ k V{4>) = ^m 2 2 (r fc ). 



(B25) 



There is no closed form solution to this. Suppose we take the standard slow-roll-like approximation except with 
H 2 w 27/3 instead of H 2 « 7/3. We find 



2 Afi 



3x/AT 



We would then write Eq. (|B23() as 



This can then solved for At. as 



where 



2 1 



3A fc k 2 /a 2 (r k )+m 



l + \ k /22k 2 /a 2 (r k )+m 

i + yi + 72/ 2 

18/ 2 

k 2 /a 2 (r k ) + m 2 x 
2k 2 /a 2 (r k ) + m 2 x ' 



(B26) 



(B27) 



(B28) 



(B29) 



28 



Now a(rfe) still needs to be determined. Since / is a function bounded by 1/2 and 1 (with the long- wavelength being 
1), we have 

0-53 < A fc < 1.2 (B30) 

where the long-wavelength limit corresponding to Xk ~ 0.53. For mx 3> m<f, ~ H, we will have mostly nonrelativistic 
contributing to the particle production. For [k/a(rk)] 2 — m x , we have Xk ~ 0.84. Hence, a more realistic range for 
X k is 

0-53 < A fc < 0.84. (B31) 

From now on, we will take this to be the uncertainty range of this parameter. 
This means that a(rk) is bounded as well. Since by Eq. I|B26J) . we have 

0.89M P < 4> < l.IMp (B32) 

where the long- wavelength limit corresponds to l.IMp . Since the Friedmann equations give 



we can write 



A ° - ^ + ± (B34) 



a(r kl ) M P V 6 y Xk 

where Aa/a = [a(rk 2 ) — afrk^/a and A<f>/Mp = [4>( r k 2 ) ~ 4>{ r kx)]/ 'Mp are both assumed to be small quantities and 
Afej w 0.53 while Afc 2 w 0.84. This results in 

A " 0.21. (B35) 



o(r fel 



To find density of particles produced along with the uncertainty, we can integrate Eq. I|B24J) assuming both A = Xk { 
and a(r) = a(rk { ) are constants independent of k and then evaluating the integral for two extreme values of Afe 4 . The 
integral is 



dkk 2 2 _ / oCrfcJ A 3 (m m x ) 3/2 /¥ /_4mx / 2 + X ki \ 

a*(2^) m ~\a(t)J 1287r2[(2 + A fcl )/3A fc J3/4V 2 ex P| m W 3Afci )■ ^ 



Hence, we can bound the number density to be 

3 



0.7 X 10^ 3 ( ^A) ( m<pmx f^ e ^-0m x /m, < n < 1.6 X 1Q-3 ( ( TO0TOx )3/2 e -4.2m x / m ^ (337) 

V a w J V a w / 

Expectedly, note that the mx dependence for mx — > of the density here differs from that of 11]. The reason 
why this is expected is that Eq. I|B9() no longer applies in that limit. We use Eq. (|B37|) in the main body of the paper 
(e.g., in Eq. 
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